Semiclassics around a phase space caustic: an illustration using the Nelson 

Hamiltonian 
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The semiclassical formula for the coherent-state propagator is written in terms of complex classical 
trajectories of an equivalent classical system. Depending on the parameters involved, more than 
one trajectory may contribute to the calculation. Eventually, however, two contributing trajectories 
coalesce, characterizing what is called phase space caustic. In this case, the usual semiclassical 
formula for the propagator diverges, so that a uniform approximation is required to avoid this 
singularity. In this paper, we present a non-trivial numerical application illustrating this scenario, 
showing the accuracy of the uniform formula that we have previously derived. 



Quantum propagators are the fundamental ingredients 
in any dynamical description of the quantum theory. 
They also provide an important tool for the study of 
the quantum-classical connection since their semiclassi- 
cal approximations can be intuitively interpreted in terms 
of classical trajectories. Naturally, therefore, they have 
been widely used in the context of the semiclassical the- 
ory. Concerned with the Correspondence Principle, Van 
Vleck [l| inaugurated this kind of study by performing 
semiclassical approximations on the quantum propagator 
in the coordinate representation (x"\exp(—iHT/H)\x'). 
According to his calculation, it can be written as a func- 
tion of classical trajectories of the Hamiltonian Hw, the 
Weyl symbol of H, connecting the initial position x' to 
the final x" , after a time interval T . Around forty years 
later, this result was improved by Gutz wilier and used 
to derive his famous Trace Formula 0, [|[ , which deter- 
mines how periodic (and isolated) orbits of Hw can be 
used to approximate the density of states of H. Notice 
that Gutzwiller's work goes beyond a semiclassical de- 
scription of dynamics. Actually, it filled a gap on semi- 
classical quantization methods since his formula applies 
to non-integrable systems, contrary to the earlier Bohr- 
Sommerfeld Q and Einstein-Brillouin-Keller [B-0] quan- 
tization rules. 

Although these seminal developments on the semiclas- 
sical theory have essentially involved approximations in 
the coordinate representation, similar results can also be 
obtained by working with other representations. In par- 
ticular, as classical states are usually points in phase- 
space, it is claimed to be natural to do semiclassi- 
cal physics using the most localized quantum states in 
phase-space, requirement accomplished by the coherent 
states [H, Q . This natural predisposition of these states 
is corroborated by the great interest in time-evolution 
problems based on phase-space pictures, which can be ap- 
preciated, for instance, in Heller's papers Moreover, 
the use of this representation has another advantage to be 
involved in semiclassical physics, namely, it can be eas- 
ily extended in order to include spin degrees of freedom 
by means of spin coherent states d, Q . This additional 
advantage can be identified already in the first paper 
that presents a semiclassical formula for the coherent- 
state propagator [ll[, since it considers both canonical 



and spin coherent states. At last, we recall that a deriva- 
tion of the Gutzwiller Trace Formula using canonical co- 
herent states can be found in Ref. [l2j|, while the analog 
for spin coherent states can be found in Refs. (l3l. Il4j]. 

In this paper, we shall focus on semiclassical approxi- 
mations of the two-dimensional coherent-state propaga- 
tor 



K(z",z',T) = (z"\exp(-iHT/h)\z') 



(1) 



The states |z) = \z x ) g) \z y ) are the coherent states that 
can be associated to a mass m subjected to a harmonic 
potential with frequencies uj r = h/(mb^.), with r assum- 
ing x or y. They are the eigenstates of the annihilation 
operator a r , namely, d r \z r ) = z r \z r ), where 



1 / q r .p r 
a r = —=[—+ 1 — 

V2 \o r c r 



and 



1 / q r | .p r 

\/2 \b r Cr 



(2) 



Here, c r = h/b r , q r and p r are, respectively, the po- 
sition and momentum operators, q r = (z r \q r \z r ), and 
p r = (z r \p r \z r ). The numbers b r and c r can be also iden- 
tified as the widths of |z) in position and momentum, 
respectively. 

According to Refs. [ll| in the semiclassical limit, 
K(z" ,z' ,T) can be written in terms of functions de- 
pending only on trajectories of the classical Hamiltonian 
H(v, u), which is achieved by calculating (z\H\z), fol- 
lowed by the replacement of z and z* by u and v, re- 
spectively. The relation between the usual classical vari- 
ables, q = (q x ,q y ) and p = (p Xl p y ), and the convenient 
variables, u = (u xi u y ) and v = (v x ,v y ), is given by 



V2 \ b r 



■ P, 



and v r 



1 / q r 
V2 \b r 



■ Pr 



Hamilton's equations for u and v are 



i dH 
h dv r 



and 



i dH 
h du r 



(3) 



(4) 



The trajectories involved in the semiclassical evaluation 
of Eq. ([TJ must obey the boundary conditions 



u(0) = z' and v(T) = z 1 ' 



(5) 
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which imply that q and p are complex, in general. Oth- 
erwise, both initial and final phase-space points would 
be fixed by the input, so that, generically, there would 
be no trajectory satisfying so many restrictions. This 
is the reason for the change (z*,z) — > (v, u). Once we 
have found such a contributing trajectory, we evaluate 
its complex action 



S(z"*,z',T) = 



— (u • v — u • v) — H 

2 y ' 



dt - A, 



where A = f [u(0) • v(0) + u(T) • v(T)], and 



d 2 H 



d 2 H 



io \du x dv x duydvy 
The scmiclassical propagator is then given by 



dt. 



K^(z"* 1 z / ,T)=Af 

traj. 



det [£S t 



(6) 
(7) 

(8) 



where Af = e 2l z l 2 l\ z I" and 



d 2 S 



dz' x dz'J* dz'^dz'^ 

d 2 S d 2 S 
dz'd 



(9) 



dz'dz': 



For non-integrable applications, it is convenient to write 
the prefactor V of Eq. (J5J in terms of elements of M, 
which is the stability matrix of the contributing trajec- 
tory: V = y/det[(i/h)S uv ] = ^/1/dct M vv , where 



8u(T) \ 
Sv(T) ) 



M uu M uv 
M V11 M vv 



Su(0) 
<5v(0) 



(10) 



Equation ((5J) is deduced by means of a quadratic ap- 
proximation around critical paths (the complex classical 
trajectories) of K(z" ,z' ,T), written in the path integral 
formalism [la, [HI (see also Refs. [U [ItJ for the one- 
dimensional case). This is the reason why we insert the 
index ^ in the symbol K. In addition, it is explicitly 
indicated by the sum in Eq. (|8]) that, in principle, we 
should consider contributions of all trajectories satisfying 
boundary conditions ([5]). At last, as trajectories depend 
just on z"* instead of z", the label z" of K is replaced 



by z" 



Some trajectories that obey Eqs. (|5|), when used to 
calculate Eq. ([8]), give origin to non-physical results as, 
for instance, probabilities greater than one. This kind of 
problem has been reported in several papers [H, EH-Hl] , 
and it is assumed that these trajectories refer to critical 
points of the path integral quantum propagator impos- 
sible to be included in any allowed deformation of the 
original contour of integration. Usually, these trajecto- 
ries are simply excluded from the calculation. A simple 
and useful rule to identify such spurious trajectories con- 
sists in writing their contributions to as e lF l h ", so 
that we can select the ones whose imaginary part of Fq, 
defined as the zero-order term of F in its fi-expansion, 



is non-negative. Otherwise, in the formal scmiclassical 
limit h — > 0, their contributions (~ e -i™\F \/h^ ^q^^ 
produce a non-physical |i^ 2 '|. Notice that, although S 
may have terms on ft, it is a good estimate to think 
of Im[Fo] = Im[<S] — ft In AT = (for a careful discus- 
sion about the ft-depcndcncc of each term of Eq. ([5]) , see 
Ref. [HI). 

It should be mentioned that the abrupt removal of 
a contribution from Eq. (JHJ) is a manifestation of the 
well-known Stokes Phenomenon [231425T ] . Generically, 
it appears when an analytic function (Eq. (TTJ), in our 
case) is asymptotically (ft — > 0) approximated by a 
multi-valued function [the sum of exponential contribu- 
tions ©]. Stokes Phenomenon refers to the fact that 
the proper choice of a branch in the approximating func- 
tion is domain-dependent. We emphasize that the sudden 
change in the form of the approximating function does 
not represent its (numerical) discontinuity. Actually, it is 
needed in order to assure the continuity manifested in the 
function represented. The criterion concerning the sign 
of J-q presented earlier combined with considerations on 
continuity shall be, therefore, our basis to decide if a 
trajectory should contribute to the propagator or not. 
Finally, as a generic manifestation of asymptotic approx- 
imations, the phenomenon is quite often in semiclassical 
physics. Apart from the cases cited earlier concerning 
the coherent-state propagator, it can be observed, for 
instance, in the WKB method 0, [2fj| and also in the 
propagator in the momentum representation p7j . 

Besides the problem of spurious contributions in 
Eq. ©, it may appear trajectories for the which the pref- 
actor V diverges. The point where it happens is called 
phase space caustic (PSC), and it is caused by the co- 
alescence of contributing trajectories. From the mathe- 
matical point of view, it arises because second order cor- 
rections of the expansion of K(z", z', T) around the clas- 
sical trajectory vanish. To avoid this problem, we need 
to develop improved approximations where further cor- 
rections are considered. As it arises from the approach 
performed and not because of the trajectory itself, we 
point out that, in this case, there is no reason to exclude 
a trajectory from the calculation. We shall return to the 
treatment of this issue opportunely. 

In order to illustrate this rich scenario of spurious tra- 
jectories and PSC's, we apply Eq. © to 



H 



1 



(Pt+fy) + [% 



Qx 
2 



A* -2 



(11) 



known as Nelson Hamiltonian, which has been studied in 
both Classical 2^, 29[ and Quantum Mechanics [10, Hl| . 
Actually, in Ref. [15(, we had already used this system 
to study the applicability of Eq. ©. Now, we revisit 
this work in order to deal with the problem of PSC 
there presented, but not solved. As well as we previously 
did [15[ , we shall restrict the application to the case where 
z' = z" = z. In addition, we shall also use b x = b y = 0.2, 
/i = 0.1 and h = 0.05. By doing so, five numbers (q X7 q y , 
p x , p y , and T) become the input parameters to calculate 
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FIG. 1. Panels (a)-(c) show the contour plots of |i0 2 '| in the (T,q x ) plane. The individual contribution of family /„ is shown 
in panel (a), while the one of family /j is shown in panel (b). The inset of panel (a) shows the contour curve \K {2) \ = 0.33 for 
families f a (solid line) and (dashed line) superimposed. Panel (c) presents IJf' 2 '! evaluated with both families. Panels (d) 
and (e) show the contour plot of To for families f a and fb, respectively. While f a has no point where To < 0, for /(, practically 
the whole region inside the contour curve To = 10 -3 has To < 0. Panel (f) combines the plots (a) and (b). In all contour plots, 
including those of Fig. [21 some contours were highlighted to facilitate a comparison among them. For panels (a)-(c) and (f), 
the difference between two subsequent curves is 0.03 [the same for Fig. [2ja)-Fig. [2jd)] . Axes are the same for all graphs. 



A'( 2 ). We then define p x = |p| cos 9, p y = |p| sin#, and 




holding E = 0.5 and 9 = 140°. Thus, for a given input 
pair (T, q x ), we select q y by the rule q y = 2q x /3, so that 
the last undetermined parameter p is solved by the last 
equation. We point out that this set of parameters was 
chosen in order to find a region containing a PSC, but de- 
scribed by a reduced number of variables, namely, T and 
q x . In the following, we show the evaluation of Eq. J8]) 
in the plane (T,q x ) for the interval 0.2 < q x < 1.0 and 
7.0 < T < 8.0. 

For all points (T,q x ) considered, we found two con- 
tributing trajectories to (32|. Based on continuity 
criteria, we can distinguish two families of such trajecto- 
ries, f a and fb- The individual contribution of each fam- 
ily to Eq. © is shown in Fig. QTa), f° r fa, and Fig.QTb), 
for /(,. Although a clear vertical cut line appears in the 
plots, we point out that the combination of both families 
gives origin to a continuous two-branch surface exhibit- 
ing a Mobius strip structure. In the inset of Fig. HJa), 
we demonstrate this property: By circulating the con- 
tour curve \K^ \ = 0.33 of f a , to avoid the discontinuity 
at the cut line, one should change to the contour curve 
|if( 2 )| = 0.33 of fb- Then, if one continues to follow this 



curve, one arrives again at the cut line, where one can 
continuously return to family f a , closing a cycle of two 
turns. In Fig. QJc), wc evaluate |_KT( 2 )| including both 
families. Notice that these three plots present a sharp 
peak that demands investigation. In Fig.[5{f), therefore, 
we plot the results of Figs.[TJa) andQJb) just for the line 
q x = 0.58, where we clearly identify a divergent behavior. 
As shown in the inset of this figure, it appears because 
of the presence of a PSC, point where V~ Y goes to zero. 

Apart from the region under influence of the PSC, 
which can not be properly evaluated by Eq. ([8]) , the ques- 
tion that naturally arises is about which plot satisfacto- 
rily approaches the equivalent full quantum mechanical 
calculation. In order to answer this point, we plot in 
Fig. QTd) and Fig. [lie) the value of for f a and 
respectively. From them, and according to the criterion 
concerning the sign of Tq defined above, we conclude 
that there is no reason why to exclude f a from Eq. ([S]). 
On the other hand, for fb, there is a large region where 
Tq < 0, implying that this family can not be used to 
evaluate |if^ 2 )| in this region. Then, if we assume that 
only f a contributes to this particular region and impose 
continuity in the whole plane (T,q x ), we find the result 
shown in Fig. [Tf f ) . In practice, to plot this graph, ex- 
cept for the region where we know that fb should not 
be included, wc span the whole plane (T, q x ), comparing 
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the results of Figs, QJa) and [1(c) and selecting the one 
which optimizes continuity. Figure QJf), in this sense, 
is the better we can do by using Eq. dHJ. A close look 
at Figs. [1(a), [T(c), [He), andQJf), however, reveals that 
family ft, was excluded from the calculation even where 
J-'o is non-negative. It could be seem as an illegitimate 
procedure, but, as discussed earlier in the present paper, 
supported by the Stokes Phenomenon, we are allowed to 
do this. 

We now face the problem of the PSC. There is no 
other solution to this issue unless to revisit and to im- 
prove the approximations performed in Eq. (TTJ). Two 
ingredients are crucial to accomplish this task: Maslov's 
method [33l - l35| and uniform approximations [23|, l36T - [40j . 
Generically, the first one consists in working with two 
conjugate representations of the same semiclassical ob- 
ject, so that if there is a singularity in a given representa- 
tion one changes to other. Then, by transforming back to 
the original representation including further corrections, 
the singular point can be circumvented. The method has 
the advantage of avoiding the extremely complicated cal- 
culation originated by the direct implementation of third 
order corrections in the path integral representation of 
Eq. (JTJ) . Uniform approximations, in this case, are use- 
ful to perform the integral involved in the last step of 
the Maslov method. They enable us to map the compli- 
cated integrand into a simpler one, but having a similar 
structure of saddle points. Essentially, as in the asymp- 
totic limit h — > saddle points concentrate useful infor- 
mation about the integral, the mapping assures a good 
accuracy. We recall that the combination of both meth- 
ods were already used to deal with problems similar to 
PSC's, namely, to deal with the break-down of quadratic 
approximations. It was used, for instance, to avoid the 
turning point divergence in the WKB method (2(|[34| and 
also to treat caustics in the coordinate propagator [4l| . 
In Rcf. j4l| . in particular, Schomerus and Sieber derived 
extensions of the Gutzwiller formula for the case of co- 
alescent orbits. We point out that, while in Ref. [4l[ 
the methods were applied to perform the trace and the 
Fourier transform of the coordinate propagator (opera- 
tions which transform the propagator into the uniform 
density of states), our approximation is exclusively per- 
formed in order to find a uniform formula for the propa- 
gator itself. 

Following these ideas, we studied a conjugate represen- 
tation for the coherent states [42jj. From it the Maslov 
method can be applied so that a uniform approxima- 
tion for K, valid for regions close and far from PSC's, 
can be easily obtained. This task was already performed 
for both ID [43| and 2D coherent states [iij . with some 
simple applications presented in Ref. |45| . The prescrip- 
tion that we achieved to properly evaluate the propagator 
around a PSC includes, firstly, finding the two contribut- 
ing trajectories to K^ 2 \ Then, from their actions iS a and 
Sb, we calculate 

A=^(S a +S„) and B = [^{S b - S a )] V \ (13) 



which can be directly used in the uniform formula [45j | 

if (uil) (z"*,z',T) = iv^c^-CB) + (B)] e A , (14) 
where c\ = (h\, — h a )/\/B and = h a + hi,, with 

h a ,b = VrVBAdetM™)!.,,, e^ b . (15) 
The function f, (£) is the well-known Airy's function, 

m = ^l c ^{i^t + ^j) d t, as) 

where the index j refers to three possible paths of inte- 
gration Cj , related to three different Airy's functions |4(J ■ 

When used for the present application, the three pos- 
sible solutions for the uniform formula (|14p can be easily 
organized as three continuous and distinct solutions in 
the plane (T,q x ), as shown in Figs. [Ha), [2(b) and [2(c). 
Comparing them with the results of Fig. [TJ we realize 
that Fig. [2(a), Fig. [2(b), and Fig. [2(c) refer, respec- 
tively, to the uniformization of Fig. [1(a) , Fig. [1(b) , and 
Fig. [1(c). Since Fig. [2(a) agrees with Fig. [1(f) far from 
the PSC, we elect it as the final semiclassical result. In 
Fig. [2(d), we show the full quantum mechanical result, 
and, in Fig. [2(e), the relative error | \K\- |A'( un > | \/\K\ be- 
tween Fig. [2(a) and Fig. [2(d). Notice that the adopted 
approximation satisfactorily agrees with the exact result: 
For almost all points, the error is less than 5%. At last, 
for an additional comparison, in Fig. [2(f), wc plot in the 
same graph the two individual contributions for K^ 2 \ 
the exact, and the uniform result of Fig. [T(a), just for 
the line q x = 0.58. We emphasize that the choice of the 
branch of Eq. (TH)) can also be seem as a manifestation 
of the Stokes Phenomenon. Notice also that, once the 
singularity has been removed, continuity is more evident 
and the choice becomes easier. 

In this article, we briefly reported a rich numerical 
study on semiclassical approximations of the coherent- 
state propagator. Analogously to the most famous semi- 
classical approximation based on second order expansion, 
namely, the WKB formula, also suffers the problem 
of non-physical solutions and singularities in its prefac- 
tor. The first problem can be handled by excluding spuri- 
ous trajectories using continuity (physical) criteria. This 
elimination, however, does not solve the problem of PSC 
(equivalent to the turning point divergence in the WKB 
method). Combining the Maslov method with a conju- 
gate for the coherent-state representation, we developed 
a uniform formula for K, which has shown to agree with 
the full quantum mechanical calculation. Curiously, in 
the uniform approach, spurious trajectories of be- 
come crucial to calculate K^ un ^ . Concerning the decision 
about which solution of if( un ) should be chosen, it can 
be done by imposing continuity and the fact that K^ un ^ 
should agree with in regions far from PSC's. Fi- 

nally, we emphasize that the numerical example studied 
in the present paper figures among the worst scenarios to 
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FIG. 2. Panels (a)-(c) show three continuous and distinct contour plots for the solutions of |i£'' un - ) | in the (T,q x ) plane. 
Panel (d) shows the contour plot of the exact \K\, and panel (e) the relative error between panels (a) and (d). In panel (f), we 
present, just for q x = 0.58, the results of Fig.[2{a) (gray solid line), Fig.[2{d) (black solid line), Fig. [TJa) (dash-dotted line), and 
Fig- Dlk) (dashed line); Its inset shows IT 3 ! -1 for f a (dash-dotted line) and /(, (dashed line), for the same points. For panel (e), 
the difference between two subsequent curves is 0.01. Axes are the same for all graphs, except for panel (f) where they are 
explicitly shown. 



evaluate semiclassical formulas, since we approach very 
close to a caustic. In spite of the adverse conditions, the 
systematic use of techniques concerned with both stan- 
dard second order approximation and uniform approxi- 
mation has shown to be quite satisfactory to deal with 



this situation. 
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